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Abstract 

In this paper, we study the matrix model proposed by Berenstein, Maldacena, 
and Nastase to describe M-theory on the maximally supersymmetric pp-wave. 
We show that the model may be derived directly as a discretized theory of su- 
permembranes in the pp-wave background, or alternatively, from the dynamics 
of DO-branes in type IIA string theory. We consider expanding the model about 
each of its classical supersymmetric vacua and note that for large values of the 
mass parameter fi, interaction terms are suppressed by powers of fi~^, so that the 
model may be studied in perturbation theory. We compute the exact spectrum 
about each of the vacua in the large /i limit and find the complete (infinite) set 
of BPS states, which includes states preserving 2, 4, 6, 8, or 16 supercharges. 
Through explicit perturbative calculations, we then determine the effective cou- 
pling that controls the perturbation expansion for large /x and estimate the range 
of parameters and energies for which perturbation theory is valid. 
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1 Introduction 



At the present time, the most promising candidate for a microscopic definition of eleven 
dimensional M-theory is the large limit of the quantum mechanics obtained from 
dimensionally reducing ten dimensional U{N) SYM theory to 0+1 dimensions, with 
Hamiltonian 

Ho^RTr Qn^-l[X^,Xs]^- i^Vf^A,*]) A,B^l,---,9. (1.1) 

The initial evidence for this suggestion came from the idea that the quantum com- 
pletion of eleven dimensional supergravity should be realized as a theory of quantum 
supermembranes just as supergravity in ten dimensions is given its quantum completion 
via superstrings. Unlike the superstring case, the worldvolume supermembrane theory 
is an interacting theory and the theory contains no tunable coupling constant; however, 
de Wit, Hoppe, and Nicolai [1] showed that in the light cone gauge, the supermembrane 
action can be regularized to give the quantum mechanics (1.1), providing a possible 
avenue for directly studying the quantum supermembrane. It was soon realized that 
in contrast to the discrete spectrum of states obtained for a quantum superstring, the 
spectrum of the matrix model (1.1) is continuous, and at the time this appeared to be 
a rather troubling feature of the quantum supermembrane. 

The matrix model was given new life following the realization that a quantum com- 
pletion of eleven-dimensional supergravity should arise directly from string theory as 
the strong couphng hmit of type IIA strings. Based on this connection. Banks, Fischler, 
Shenker, and Susskind (BFSS) [2] argued that M-theory in the "infinite momentum 
frame" should be described by the limit of a large number of DO-branes in type IIA 
string theory at low energies and weak coupling. Remarkably, this system is described 
again by the quantum mechanics (1.1), so essentially the same conclusion is reached 
from a very different approach. Happily, the DO-brane point of view provided a very 
natural explanation for the continuous spectrum of the matrix model: the theory should 
be interpreted as a second quantized theory containing both single and multiparticle 
states, rather than a first quantized theory. This point of view also provided a precise 
interpretation for the finite N matrix model, as the DLCQ of M-theory with A^ units 
of momentum in the compact direction [3]. 

Since the work of BFSS, matrix theory has passed numerous tests (for a recent 
review, see [4]), however an honest quantum mechanical study of the model seems 
extremely difficult, due to the fiat directions in the potential (resulting in the contin- 
uous spectrum mentioned earlier), the difficulty of distinguishing single-particle and 
multi-particle states, the lack of a tunable coupling constant (interactions are of order 
one), and the large number of degrees of freedom in the large N limit. Indeed, even 
determining the normalizable ground state wavefunction for the case N — 2 (proven 
to exist in [5]) is an unsolved and notoriously difficult problem. 

Recently [6], it has been shown by Berenstein, Maldacena, and Nastase (BMN) that 
the matrix model (1.1) is actually part of a one-parameter family of matrix models with 
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32 supercharges given by 



+ i^e^^'X,X,X,) . (1.2) 

The matrix model for general /i was proposed by BMN to give a description of the 
DLCQ of M-theory on the maximally supersymmetric pp-wave background [7, 8, 9, 10] 
of eleven-dimensional supergravity^ 

9 / 3 2 9 2 \ 

els'" = -2dx+dx- + y dx^dx^ - y —x'x' + y ^x^x" ] dx+dx+ 

i^l23+ = 1^ (1.3) 

In the large N limit with fixed = ^ and fixed /x, this matrix model should then 
describe M-theory on the pp-wave background in a sector with boost-invariant momen- 
tum parameter fip^ . 

As pointed out in [6] , the pp-wave matrix model for non-zero avoids a number of 
difficulties associated with the // = case. Firstly, the quadratic terms in the potential 
remove all fiat directions, and yield a potential that becomes infinite in all directions. 
This implies a discrete spectrum for the massive matrix model and, as shown in [6], 
an isolated set of classical supersymmetric vacua. These vacua may be associated with 
all possible ways of dividing up the momentum into different numbers of membranes, 
from a single membrane carrying N units of momenta to N membranes each carrying a 
single unit of momentum. Since these vacua are separated by potential barriers, we see 
that the massive matrix model (at low energies) also has a natural distinction between 
single and multi-membrane states. 

In this paper, we study the massive matrix model and explore the consequences of 
these and a further simplifying feature: that the existence of a tunable parameter in 
the matrix model opens up the possibility of a weakly coupled regime which can be 
studied in perturbation theory. 

We now provide a concise summary of our results as we outline the remainder of 
the paper. 

Before proceeding to study the model we first clarify its origins.^ In section 2 we 
show that the massive matrix model arises directly as the regularized theory of light- 
cone gauge supermembranes in the maximally supersymmetric pp-wave background. 
In section 3, we note that the model may also be derived from the DO-brane perspective, 
although here one has to assume that the weak background DO-brane action continues 
to be valid even in a region of high curvature, perhaps a consequence of supersymmetry. 

^We have used conventions in which = (t ± x)/y/2, different from those in BMN. 

^The derivation provided in [6] was somewhat indirect, showing that the = 1 model matches 
the action for a superparticle in the pp-wave background and then writing an extension of the theory 
to general N consistent with supersymmetry. 
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In section 4, we proceed with a direct study of the massive matrix model. We begin 
by reviewing the action and its symmetries and identify exphcit matrix theory ex- 
pressions for the generators of the superalgebra, verifying that the (anti)commutation 
relations are satisfied correctly. We note that the appearance of rotation generators 
along with the Hamiltonian in the anticommutator of supersymmetry generators allows 
BPS states carrying angular momentum.'^ After these generalities, we review the classi- 
cal vacua of the theory, collections of fuzzy spheres corresponding to sets of membrane 
giant gravitons [11, 12] carrying various fractions of the total light-cone momentum, 

. We expand the matrix model about each of these classical vacua and note that 
in each case we obtain a quadratic theory with interactions suppressed by a coupling 
We conclude that for fixed values of the other parameters, the model may be 
studied in perturbation theory for sufficiently large 

In section 5, we consider the theory in the // — > oo. Here, the different vacua become 
superselection sectors separated by an infinite potential barrier and interactions vanish. 
For each vacuum, we explicitly diagonalize the quadratic action, obtaining the precise 
oscillator masses and thus the exact spectrum of the theory. We find that the oscillator 
masses are proportional to jJL and completely independent of N . In fact, N plays the role 
of an ultraviolet cutoff leaving a finite tower of oscillator masses that becomes an infinite 
tower for N ^ oo. This fits very well with the regularized supermembranc picture, 
and in particular, we may interpret the large N limit of the spectrum obtained from 
the single fuzzy sphere vacuum as the quadratic spectrum of a single supermembranc 
in a large /ip'^ limit of the pp wave background. From the exact spectrum, we are 
able to exphcitly identify all BPS states of the theory and find infinite towers of states 
preserving 2, 4, 6, and 8 supercharges in addition to the 1/2 BPS vacuum states. 

In section 6, we move away from the /i = oo limit and investigate perturbation 
theory about the quadratic actions. We expand the interaction Hamiltonian in terms 
of the canonically normalized oscillators, and find that the parameter suppressing the 
interaction terms is different for different vacua. In particular, the coupling for the 
A'^- membrane vacuum remains g = {R/ n)^ while the coupling for the single-membrane 
vacuum is reduced to g = (R/ {fiN))^ . This latter quantity, equals to r~^/^ where r is 
the radius of the membrane sphere, remains finite in the M- theory limit of fixed fi and 
fixed = N/ R, giving the coupling constant for the continuum field theory living on 
the single spherical supermembrane. 

Since it is often the case that bare couplings become modified to some effective 
coupling in theories with many degrees of freedom, we then compute explicitly the en- 
ergy shifts for certain non-BPS excited states of the single-membrane and A^-membrane 
vacua. In both cases, we find that the leading corrections to the energy are suppressed 
by an effective coupling \ = g^K relative to the energy in the free theory, where g 
is the bare coupling for the corresponding vacuum and K is the multiplicity of the 
representation (number of coincident branes). Thus, the leading perturbative correc- 
tion in the irreducible vacuum is small when g"^ = l/{fip^)^ ^ 1 while a small leading 
correction for the X = vacuum requires a much stronger condition, N'^/^np'^Y 1- 

^An example of such a state has been given in [24]. 
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Figure 1: Regime of validity for perturbation theory (conservative) 



For fixed values of the parameters, we then estimate the energy scale at which 
perturbation theory will break down (i.e. the interaction Hamiltonian becomes of 
the same order of magnitude as the quadratic piece). For generic directions away 
from the singlc-mcmbranc vacuum, we find that interactions become important at an 
energy E/ ^ ~ = (nN/R)'^ but there are certain directions (corresponding to 

the large membrane shrinking and a small membrane growing at the origin) for which 
interactions become important at E/^ ~ l/{g'^N) = j^i^jiN/RY. We expect that a 
valid perturbative regime will certainly exist when this energy is significantly greater 
than the energies of the lowest excited states, g^N <C 1, as depicted in figure 1. On the 
other hand, our direct perturbative calculations indicate small perturbative corrections 
for low energy states under the much weaker condition g"^ <^1. 

In section 7, we discuss the implications of our results for the large N limit of the 
model. Our calculation of the leading energy correction indicates that perturbation 
theory near the single membrane vacuum may be valid even in the large limit as long 
as /xj?"*" 3> 1. On the other hand, the more conservative estimate based on the height 
of the barrier between vacua suggests that we always exit the perturbative regime in 
the large N limit defining M- theory (as in figure 1). In any case, our results explicitly 
demonstrate that in this limit, the theory develops finite zero-energy fiat directions 
corresponding to decreasing the radius of any membrane while increasing the radius 
of any other membrane (or nucleating a new spherical membrane at the origin). This 
indicates the existence of massless singular BPS domain walls, or tensionless strings, 
in the supermembrane world-volume field theory. 

Finally, we offer some concluding remarks in section 8 and technical results in a set 
of appendices. 

2 The matrix model from supermembranes 

It is well known that the usual Matrix theory action can derived as a discretized version 
of the light-cone gauge supermembrane action in flat space, where Poisson brackets on 
the membrane worldvolume are replaced with matrix commutators [1] . In other words. 
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the invariance of the supermembrane action under area preserving difFeomorphisms is 
translated to the U{N), N ^ oo gauge invariance of the dimensionaUy reduced SYM 
theory which defines the matrix model. 

In this section wc show that the same story applies in the case of the pp-wave back- 
ground. That is, a discretized version of the light-cone supermembrane action on the 
maximally supersymmctric eleven dimensional pp-wave background exactly reproduces 
the matrix model proposed by BMN. 

To start with, we borrow the supermembrane action from [13]: 

S[Z{a)] = / d'a ^J-G{Z{a)) + / C , (2.1) 

where Z^{a) = X^{a),9°'{a) are the superspace embedding coordinates and a\i = 
0, 1, 2 are the world-volume membrane coordinates. Here, 

G{Z{a)) = det(nfn^^G^.) , (2.2) 

where 11^ — ^^E^ is the pull back of the supervielbein to the membrane world- 
volume and Gij,i, is the background metric, which in our case is the pp-wave metric 
(1.3). The three form C field is 

C = ^e" A A e* Crst + Cp , (2.3) 

where = e^AX^^ and el are the vielbeins. The Cp term, which denotes the fermionic 
contributions to the C superfield, will be presented explicitly in section 2.2. 

We begin with the bosonic terms and then consider the fermionic terms in the 
supermembrane action. 



2.1 Bosonic terms 

In order to obtain the matrix model from the membrane action one should fix the light- 
cone gauge. For that it will be convenient to use the coordinates {X~^ , X~ , X"^) with 
A — 1, ...,9. The light-cone condition is then X~^ = r, while the other worldvolume 
coordinates are denoted by cr'', r = 1,2. One should note that to fix the light-cone 
gauge consistently we should also impose the following constraint on the phase-space 
variables 

PAdrX^ + P_drX- = , (2.4) 

where Pa,P- are the momentum conjugate to X"^ and X~, respectively. 

Upon fixing the light-cone gauge, the bosonic part of the total Hamiltonian becomes 

[13] 



Hbos — 



(fa 



-Ca- \^GA^f + \{e-drX^dsX^r 



C I ^ 



e'''drX^dsX^PAC+^B + C-C+J\) (2.5) 
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where 



Ca = -e^'drX-dsX^C.AB + i^e'-'drX^dsX^CABC 
C± = -e^'' drX^dgX^ C±AB 
C+_ = e''drX-dsX^C+^A . 
Now let us concentrate on the pp-wave background, where 

G — = Ga- = 0, G+- — —1, Gab — ^ab 



(2.6) 



(2.7) 



and 



C_ = C+_ = = 0, C+ = ie^^a,X^a,X^C+^B = |e,,fce^^a,X^a,X^X^' . (2.8) 

In this background the hght-cone Hamiltonian above simphfies significantly and we are 
left with 



Hi 



bos 



\pl + \{x\ x^Y + \{x\ x'Y + \{x\ x^Y 



+ 1 ( [^f^x'f 



2 V 3 

where i, j = 1, 2, 3, a, 6 = 1, 2, • • • , 6, p"'" = — P_ and 



• (2.9) 



2.2 Fermionic terms 

To get the fermionic terms we have to study the supervielbein and the tensor field 
backgrounds to all orders in 9. For coset spaces (of which the pp-wave is an example), 
the exact expressions for the supervielbein E^,'^ in terms of the component worldvolume 
and background fields have been worked out in [14], and are given by^ 



16 



1 



^1 (2n+l)! 



15 



\ (2n + 2)! 



(2.10) 



(2.11) 



^Thc B = (r, a) indices arc tangent space indices and the curved indices are N = {fx, a), where the 
first entries are bosonic and the second entries are fermionic. 

^The coset space approach was also used in [15, 16] to derive the action for hght cone type IIB 



string theory in the pp-wave background. 
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In these expressions, 

e'' = e^dX^" (2.12) 
u''' = ulfdX^' (2.13) 

De = rf^ + e''Tf"W,t„„-^a;^T,,^ , (2.14) 
where lu^^ is the spin connections, Fstuv is the field strength of the three form C and 

T/*"'' = ^|^(r/*"" - 85,[T*""1) . (2.15) 
Finally, the matrix Al^ is defined as 

{M'Tb = liT^'^^^efF^tUer), - ±^{T,ser(e[r'''''""Ftu.^ + 24r,„F^^*"])^ . (2.16) 

In the pp-wave background, 

e-_ = e\ = l, e\^5\, = ^ , (2.17) 

and all the other components are zero. Then, the only non- vanishing component of the 
spin connection is 

cu;^ = \dAF^ . (2.18) 

To work out the explict form of the fermionic fields we should choose a specific 
representation for V matrices. The 32 x 32 F matrices can be written as 



''-^{a^ oj . r-V2(^ '7j , r^4; (2.19) 

where the are 16 x 16 Euclidean 5*0(9) gamma matrices. ^ We can also use 
K-symmetry to fix the (light-cone) gauge: 

F+^ = . (2.20) 

In this gauge we have the advantage of significant simplifications such as 

^F+ = ev+'de = er+^de = ev+^^oe = ev^de = ev^^ae = o, • • • (2.21) 

One can solve the gauge condition (2.20) exphcitly: 

Ml)' »=2^(-*' »)• P.22) 



e 



where ^' is a 16 dimensional (Majorana) spinor. The normalization for ^> have been 
chosen so that 



*^Note that r± = ^{V^ ± T^O), and 6 = iO^T'^. 
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To work out the action explicitly, we first note that for the pp-wave case, as discussed 
earlier, the only non-trivial background field component appearing in is F+123 = 
_^-i23 Therefore, using the gauge condition (2.20) one can show that 

= . (2.23) 

That is, the only non-vanishing fermionic contributions come from the expansion of 
the supervielbeins to the leading order in Q. 

Now let us focus on the fermionic contributions coming from the \J —G terms in the 
action. For these, we need to work out H^'s: 

n[ = EldiX"" + EldS'^ . (2.24) 

Using the Eq.(2.10), we have 

El = ei^erD^Q, 

El = QV'DJ , (2.25) 

where 

and E^ — —{QV^^a.- Putting this all together we obtain 

n- = 9,x- + a,x+(:^ + ^*V''*) + «*^5,* , 

n+ = nf = diX^ . (2.26) 

We can see that all the effects of the pp-wave background appear in the diX^ factor 
of n~. Hence, in the membrane Hamiltonian (2.5) G++ should be shifted to 

+ ^^^^^^^ , 

which leads to a "mass term" for the fermions. Note that the other fermionic term 
will not contribute to the light-cone Hamiltonian (because it involves first order 
time derivatives). 

The only remaining terms to be considered are the fermionic contributions to the 
Cterm, Cp [14]: 

Cf = - [ tdt eVrs Eit) A E''{t) A E'it) , (2.27) 
Jo 

where t is an auxiliary parameter and E{t)^ E^{t) are obtained from (2.10) by the shift 
9 — >• t9. As we discussed, using the gauge fixing condition (2.20), is zero and hence 

Cf^- t t^dt eVrsDe A (e'' + eev^DO) A (e^ + eevDO) . 

Jo 

In the gauge (2.20), OVrsDO is zero unless {rs} = {+A}, while OT^' DO is zero if r is -|-. 
Therefore the only surviving term is 



Jo 



Cf ^ - rdt eVr.De A A 



^ ^e'^''9r-^diX+dj9dkX^ 



(2.28) 



In summary, gathering together all fermionic and bosonic terms, we have shown 
that the total light-cone gauge Hamiltonian for membranes in the maximally super- 
symmetric pp-wave background, with properly fixed K-symmetry, is: 



Hi 



total 



I 



pi 



2 ^ 4^ -' 2 



■^e'^^{X\X^}X^ 



V{^^, - ^/i*^ V* \ . (2.29) 

Z 4 



The matrix model Hamiltonian proposed by BMN now follows immediately by 
applying the regularization prescription of [21, 1], in which we replace functions on 
the membrane worldvolume by hermitian N x N matrices, integrals with traces, and 
Poisson-brackets with commutators, i.e. 



X^{a) 
*(a) 

{, }PB 



-X' 



N 

—Tr 
R 



In the next section, we show explicitly how the matrix model may also be obtained 
from the DO-brane perspective. 



3 The matrix model from DO-brane dynamics 

In the previous section, we have seen that the matrix model proposed in [6] arises 
directly gularized version of the light-cone supermembrane action in the maxi- 

mally supersymmetric pp-wave background of eleven-dimensional supergravity, just as 
ordinary Matrix theory arises from the flat space supermembrane. 

In the usual case, an alternate approach to deriving the matrix model as a descrip- 
tion of M-theory is by showing an equivalence between the DLCQ of M-theory with N 
units of lightlike momentum [3] and a limit of type IIA string theory with N DO-branes 
in which the only remaining degrees of freedom are the lightest open strings between 
the DO-branes [2, 18, 17]. In this section, we will try to obtain the matrix model of 
BMN directly from this DO-brane perspective. 

3.1 Relation to type IIA string theory 

In the usual arguments [17, 18] , DLCQ M-theory with N units of momentum is obtained 
(or defined) by taking a compactification of M-theory on a spacelike circle of radius 
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Rs, performing a boost with parameter 



and taking the hmit i?^ ^ where the boost parameter becomes infinite leaving a 
hghthke circle of fixed radius R. The original frame, with spacelike circle of radius 
is equivalent to type IIA string theory with DO-branes in a limit Qs — ^ 0, 
a' — » in which only the low energy DO-branc degrees of freedom remain. In order 
that the energies of the corresponding DO-brane states (minus the diverging constant 
piece N/Rs) match with the fight cone energies in the DLCQ description, we must also 
rescale all dimensionful quantities in the DO-brane theory by m — > {R/Rs)m before 
taking the limit. 

In the present case, we begin with a non-flat background for DLCQ M-theory, 
namely the pp-wave metric. Applying the same arguments, we should flnd that the 
matrix model appropriate to this background is that describing DO-branes in a type 
IIA background obtained by (un) boosting and compactifying the eleven- dimensional 

pp-wave background. 

To determine the relevant llA background in the present case (where the back- 
ground flelds are not necessarily weak), note that under a boost to the speed v in a;^ 
plane, 

+ + - + ^ - 

— >\ - X^ , X — >\ - X . 

\l 1 + v ^ 1-v 

Hence, the pp-wave metric is boost covariant, and the effect of the boost is simply to 
rescale yU. In fact, up to a factor of y^, this rescaling is later removed when we rescale 
dimensionful parameters to make light-cone energies match with the energies of DO- 
brane states, so we will take both into account at once simply by replacing ji ^ ji/ \f2. 
More details of the boost and rescaling may be found in [19], where this procedure was 
carried out in detail for an arbitrary weak background. 

Taking into account the rescaling, the metric for M-theory in the "unboostcd" 
frame, which we should compactify on a spacelike circle to obtain the desired type IIA 
metric, is 

p2 p2 p2 3 9 

ds'^ = -(1 + -r)dt^ + (1 - -r)(^^a;^°)2 dtdx^° + Y,dx,dxi + ^ dXadxa , (3.2) 

^ ^ i=l a=4 

where is defined in (2.7) and the four form field is 

-^123(10) = —-^0123 = ■ (^'^^ 

Compactifying along x^^, the corresponding IIA background (in string frame) is 
dslo = goodt'^ + QAsdx^dx^ , 

-2<l> r M J 

900 ^-e 3, QAB^SABe^, A,B^l,---,9 (3.4) 
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and 



-^0123 ~ ' -^123 = ^ (3.6) 

where Co is the (time-component) of the RR-one form potential and H123 and F0123 
are the NSNS and RR-three form field strengths, respectively. 

Although the eleven dimensional pp-wave metric that we started with has zero 
scalar curvature (and its only nonvanishing component of Ricci scalar is R^^ = |/i^), 
the ten dimensional metric has a non-zero curvature. This can be understood noting 
that the effective compactification radius (e") has a non-trivial x-dependence. By 
explicit calculation one can work out the scalar curvature of the metric (3.4): 



7^=-f (1-FV4)- 



3 9 

(3.7) 

=1 a=4 



where a, b, c are positive constants. 

We can trust the above SUGRA solution only when the curvature is small, i.e. 
when <^ 4. In fact, we see that at F^ = 4, the background becomes singular, and 
is ill-defined beyond this point. Prom the eleven dimensional point of view, the reason 
for this bad behavior is that the direction becomes timelike beyond this point, so 
clearly, compactification does not make sense. Nevertheless, we may hope that in the 
region for which the type llA background is well behaved, the action describing DO- 
branes in this background may be a valid description of the regime in which the DLCQ 
M-theory degrees of freedom stay within the weakly coupled region of the pp-wave 
background. 

In this weakly curved limit, the IIA background fields can be expanded in powers 
of F^. Up to the first order in we find 

(, ^ - 

16 

^2 



4<t> 

e 3 = 
Co = 


4 
F^ 


4-F2 


900 = 


-24 
—e 3 ~ 


9ab = 


20 

6 3 dAB' 



4 

-1 - If^ 



(3.8) 



We proceed in the next subsection to determine the action for DO-branes in this back- 
ground. 



3.2 DO-brane dynamics on compact ified reduced pp- waves 

We have seen that it is only possible to relate the eleven-dimensional pp-wave metric to 
a reliable type IIA supergravity background in the region where <^ 4. In this weakly 
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curved region, we may read off the action describing DO-branes in the baclcground (3.8) 
from the results of [19], which provides the action for DO-branes at low energies and 
weak coupling in a general weak background. 

In the notation of that paper, terms arising from the eleven-dimensional metric for 
our background are 

^2 = ^^A^B^/r^ + i^^^B VC^^""^ + ^^AaBCo/"^^^) (3.9) 

where, for example ijp^^^ is the AB moment of the current coupling to (f). Using the 
expressions for the currents / given in [19] and the background fields given in (3.8), we 
obtain 

^2 = -^dAdBF\T++^^^^ + higher orders) , 

where 

R 

and the higher order terms vanish in the a' ^ limit we are taking. Plugging in 
our expression for F^, we obtain exactly the bosonic mass terms in the BMN matrix 
model. 

To obtain the terms arising from the eleven-dimensional three-form field, the re- 
quired terms in the zero-brane action are 

^3 = OaB^J^^^^ + dAC^p,lf'^^^ (3.10) 
To obtain the field strengths (3.6), we may choose potentials 

■Bij — —^ijk^kj Coij — —€ijkXk (3.11) 

Using these expressions and the expressions given for the currents in [19], we find 

S3 = uSijkiJ^''^^''^ + higher orders) 

where 

= -J-Tr {i[X\X^]X'' + -^^7^^'*^*) 
6it 8 

and again, the higher order terms vanish in the limit we are considering. Thus, the 

form-field background gives rise to both the fcrmion mass terms and the bosonic cubic 

interaction appearing in the BMN matrix model. 

We have now seen that all terms in the BMN matrix model may be reproduced by 

considering DO-branes in an appropriate weak supergravity background.^ Prom this 

point of view, it is remarkable that the model does not receive any corrections, since 

we have seen that the type IIA background becomes strongly curved and unreliable as 

we move far from the origin and in fact fails to make sense for > 4. 

''These results could also be obtained directly using the general results for matrix theory in an 
arbitrary weak background [20]. 



13 



4 Classical supersymmetric vacua 



In this section, we begin a direct study of the matrix model. We first write down 
the action in a convenient form and discuss the symmetries of the theory and their 
consequences for the spectrum. We then recall the supersymmetric vacua of the model 
and expand the matrix theory action about each of these solutions. We find that near 
each vacuum, the theory is described by a quadratic Hamiltonian with interactions 
suppressed by ^. 

4.1 Action and symmetries 

It will be most convenient to write the action in a form that makes the SO{S) x S0{6) 
symmetries manifest, so we divide the scalars into X* with i = 1,2,3 and X" with 
4 < a < 9, and write the fermions as ipja where / is a fundamental index of SU{A) ~ 
SO {6) and a is a fundamental index of SU{2) ~ SO{Z). The relation between these 
fermions and the real 16 component spinor is detailed in appendix A. 

With these conventions, the matrix model describing the DLCQ of M-theory in the 
maximally supersymmetric pp-wave background is given by^ 

+\[X\X^f + + \[X\X'^f^ (4.1) 

Here, cr* arc the usual Pauli matrices and gf j relate the antisymmetric product of two 
SU{'i) fundamentals to an 5*0(6) fundamental. 

The new terms relative to the fiat space matrix model appear in the second line. 
These include mass terms for each of the matrix theory variables, as well as a cubic 
interaction which favors matrix configurations carrying a dipole moment of membrane 
charge. 

As for the case of fiat space, the U{1) part of the theory (describing the center 
of mass degrees of freedom) decouples from the SU{N) part. In this case, the U{1) 
sector is described by a harmonic oscillator Hamiltonian with bosonic oscillators in the 
SO (6) directions of mass | and in the SO (3) directions of mass |, as well as 8 fermionic 
oscillators of mass ^. Thus, unlike the flat space case, the different polarization states 
have different masses. 

In appendix B, we review the full symmetry algebra of the matrix model and provide 
explicit expressions for the generators in terms of the matrix theory variables. The 
bosonic generators include the Hamiltonian, rotation generators in the SO (3) and 

^Throughout this paper, we set /p = 1, but we can restore Ip in any formulae using the fact that 
R and l/fj, have dimensions of length. 
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5*0 (6) directions, the central light-cone momentum P+ = ^, and generators and aa 
which are exactly the harmonic oscillator creation operators from the U{1) part of the 
theory. In addition, we have 32 fermionic symmetries. Sixteen of the generators are 
the 8 complex fermionic creation operators from the U{1) part of the theory, 

qia = :^Tr (l/ji^) 

which generate the overall polarization states. The other 16 supersymmetries are gen- 
erated by 

which include new terms linear in X relative to the flat space model. 

The supersymmetry algebra has a number of unusual features compared with the 
flat space case. First, the supercharges do not commute with the Hamiltonian but 
rather obey commutation relations 

[H,qia] = -^qia 

[H, Qia] — ^Qla 

Thus, different members of a multiplet of states generated by acting with supercharges 
will have different energies, although the energy differences will still be exactly de- 
termined by the symmetry algebra. In addition, we now have rotation generators 
appearing on the right hand side of the {Q, Q} commutation relation along with the 
Hamiltonian, 

Note that we still have {Q^'^", Qia\ oc H, so that a state is fully supersymmetric if and 
only if it has zero energy. On the other hand, unlike the ordinary matrix model, we 
now have the possibility of BPS states with non-zero energies carrying some angular 
momentum. We will see that 1/4, 3/16, 1/8, and 1/16 BPS states are possible in 
addition to the ground states which arc 1/2 BPS (preserving 16 supercharges). In this 
case, the rotation generators are not central charges since they do not commute with 
the supersymmetry generators. As a result, the representation theory of this algebra 
is quite different from the usual case, and in particular the notion of BPS multiplets is 
altered. We will discuss the BPS states in more detail in section 5.5. 

4.2 Classical supersymmetric solutions 

We now recall the classical supersymmetric vacua of the matrix model given in [6] . To 
begin, we note that the bosonic potential may be written as a sum of positive deflnite 
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terms, 



T/=-Tr 
2 



^X' + ie'^^X^X'^^'' + ]^{i[X\ X%^ + {i[X\ X']f + (X-^) 



For supersymmetric solutions, each of the four terms must vanish independently. It is 
clear that for the last term to vanish, we must take X"' = 0, while for the first term to 
vanish we must have 

where J* form a representation of the SU (2) algebra. 

Thus, the classical supersymmetric solutions are labeled by the N dimensional repre- 
sentations of SU{2), and the number is given by the number of partitions of N into 
sums of positive integers, as noted in [6]. 

Geometrically, a solution labeled by a partition {iVi, . . . , Nk} corresponds to a set 
of membrane fuzzy spheres with radii ~ The spheres are stabihzed through a 
competition between the quadratic and quartic terms in the potential, which make the 
spheres want to contract, and the cubic term which makes the spheres want to expand. 
The net effect is that any membrane in the pp-wave background expands to a preferred 
radius proportional to its light-cone momentum. Unlike the similar situation studied 
by Myers [11], the various solutions all have zero energy and are thus classically stable. 

It is particularly interesting that the classical vacuum states each correspond to a 
definite number of membranes. Thus, at low energies, one of the difficult features of 
the usual matrix model, distinguishing between single and multi- membrane states, is 
absent. 



4.3 Expansion about the classical vacua 

We now consider expanding the action about its classical vacua, that is taking 

X' ^ + Y' (4.2) 

oil 

for a given choice of J. 

It is convenient to rescale variables to remove all dependence on parameters from 
the quadratic action, so we take 

Y' J -Y\ X" J -X\ t -t (4.3) 
V /X V A* 

Note that energies will be measured in units of ii from now on; we will restore this 
factor in the important formulae. 

Plugging (4.2) into the action and performing the rescahngs, we find that the re- 
sulting action (in the Aq — gauge) takes the form 

S — S2 + S2 + S2 -|- <S'3 -|- jS'4 
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where 



St = IV {ii^H - - ^aV^/, i^p]) 



) 



(4.4) 



are quadratic actions for X, F, and and 



^3 = 




Tr 0-[J\ X''][Y\ X"] + ^[J\ Y^][Y\ Y^] - le'^^Y'Y^Y^ 

\ o o o 

Tr y^f + \[x\ x'f + x"]^) 




) 



3 



^4 = 



(4.5) 



are the interaction terms. 

We see that for large ^, the interaction terms are suppressed by powers of relative 
to the quadratic action, so we should be able to study the low energy excitations about 
these classical vacua using perturbation theory in ^. In particular in the ^ — > oo limit, 
the theory divides into superselection sectors labeled by iV-dimensional representations 
of SU (2) and each section will be described by a set of bosonic and fermionic harmonic 
oscillators. Our next step will be to solve the quadratic theory explicitly in this limit 
before considering the effects of interactions in section 6. 

5 Exact spectrum in the large ji limit 

In the limit of large ix for fixed N and i?, the various classical vacua become separated 
by infinite energy barriers and may be treated as superselection sectors, each described 
by a quadratic Hamiltonian. As a first step towards understanding the theory, we would 
now like to explicitly diagonalize the quadratic actions in order to find the complete 
spectrum in the limit of large fi. 

5.1 Expansion in terms of matrix spherical harmonics 

To proceed, we note that many of the terms in the quadratic action are written in 
terms of commutators with the J* matrices. 



It will be convenient then to expand our fluctuation matrices in a basis which behaves 
simply under this commutator action. To do this, note that the set oi N x N matrices 



[J\A] 
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forms an iV^-dimensional reducible representation of SU{2), with the SU{2) genera- 
tors acting as above. It is not hard to see that this representation decomposes into 
irreducible representations as 

Ar2 = 1 © 3 © . . . © (2Ar - 1) , 
that is, integer spins up to j = A?" — 1. We may therefore expand a,n N x N matrix as 

AT-l j 
j=0 m=-j 

where Yj^ transform in the irreducible spin j representation. For Hermitian matrices, 
we must have the reality condition 

a*^ = {-ira,_^. (5.1) 

Explicit expressions for these matrix spherical harmonics may be obtained by writ- 
ing the usual spherical harmonics as symmetric polynomials in Xi with xl + x^ + xl = 1 
and replacing x^s with Js (up to normalization). In this way, we find that the spin j 
representation may be written as products of j matrices whose indices are contracted 
with symmetric, traceless tensors. A detailed discussion may be found in [21]. 

For example 



^^°-V(7V^"^' ^^^-"V(7V^"^^'--- ^^-^^ 

where — ± iJ"^. In general, we will choose normalizations so that 

Tr {Yj,^,YjJ^N5j,j6^,m- (5.3) 

With these definitions, we may use the usual properties of angular momentum gener- 
ators. 



[J3, J = mYj^, [J+, Yj^] = \J{j - m){j + 1 + m)Yj ^+1 

[J\ [J\ yjm]] = j{j + l)Yj,n [J-, Yjm] = {j + m){j + 1 - m)Yj m-i (5.4) 

We now proceed to diagonalize each of the quadratic actions above. 
5.2 Spectrum of the irreducible vacuum 

We begin by considering the single-membrane vacuum corresponding to the N dimen- 
sional irreducible representation of SU(2). 
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Bosons in the SO (6) directions 

Prom the action above, we see that the quadratic potential for X" may be written as 



232 

In this form, it is obvious that the mass matrix is diagonal in the basis of irreducible 
representations of SU{2), 

X°- — r" Y 

and the eigenvalues are 



1\\2 

each with degeneracy 6. The quadratic Lagrangian becomes 



where j runs from to — 1 and m runs from —j to j, as usual. 
Bosons in the SO (3) directions 

For the oscillators, we may rewrite the quadratic potential as 

2 3 

Thus, if we solve the eigenvalue problem 

Y' + ie'^''[J^,Y''] = XY' 

the masses will be given as = (|A)2.^ 

To find the eigenvalues and eigenvectors, we write 

Making this substitution and writing out the three equations corresponding to i 
+, — , 3, we find a set of equations 



(A - m)y-^_^ = + 1) - m(m - l)y]^ (5.5) 

- ^)yjm = + 1) - + i)i//m+i - + 1) - - ife-i 



^Hcrc, the matrix transforms in the tensor product of the spin 1 representation of SU{2) (from 
the exphcit vector index) with the sum of the spin j representations. The eigenvectors are precisely 
the irreducible representations, corresponding to vector spherical harmonics. 



19 



Note that we have a separate set of equations for each value of m between m = —j — I 
and m — j + 1}^ Replacing and y~ in the third equation using the first two 
equations, we find an equation for the eigenvalues, 

A(A + j)(A-(j + l)) = 

so we may have A = 0, A = — j' or A = j' + 1. Using the first two equations, we can 
read off the eigenstates for these three eigenvalues. For A = 0, the eigenstates are 

Y' = [J\Yjm] l<j<N-l, -j<m<j 

for a total of N"^ — 1 zero modes. These correspond to fluctuations in the directions 
of the gauge orbit and not to physical zero-modes. To see this, note that under an 
infinitesimal gauge transformation. 



-r 



so that any fiuctuation of the form 



Y' = [J\A] 



corresponds to a gauge transformation. Since [J^Yoo] = 0, we have N'^ — 1 instead of 
A*"^ gauge orbit directions, as we found above. 

For the A = —j eigenstates, corresponding to oscillators with M = |, we find 



— Im 
(-^j—lm 



{j + m){j + 1 + m) 



i(2j + l) 



j m+1 



^j—lm 



U - 'm){j + 1 - m) 



j(2j + l) 



y 



j m-1 



U + m){j - m) 



VN \ j{2j + 1) 



Y. 



(5.6) 



where < j < N and —j<m< j and we have the reality condition a*^ — 
For the A = (j + 1) eigenstates, corresponding to oscillators with M — we find 



y+ = 
y- = 

Y' = 



Pj+i rn U - + 1 - m) 
VTV ^ (j + l)(2j + l) 

Pj+im {j + m){j + l + m 



y 



j m+1 



■y 



y/N \ (j + l)(2j + l) 

Pj+l m 0' + 1 + "00' + 1 - m) 
^\ (j + l)(2j + l) 



3 m-l 



Y 

^3 m 



(5.7) 
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For m = ±j, one of the equations vanishes while for m = ±(j + 1), two of the equations vanish. 



20 



where < j < and — j — 1 < m < j + 1 and we have the reahty condition 
In terms of a and /3, the quadratic Lagrangian becomes 

Y ^-2 1 . 2 1 J + 1 2 2 



^ 2" 2'3' 



Fermions 

The quadratic fermion potential may be written as 

V^^^^'^{\^a + la'AJ\M (5-8) 

where we have suppressed SU{4) indices. Thus, we would like to find the eigenstates 
of the equation^^ 

in terms of which the fermion masses will be 
Making the expansion 

Ya Ya jrn 

and inserting into the eigenvalue equation, we find two equations 

(A-m)V'f = ^J{j-m){j + l + m)^|;i_^^' 
(A + m + l)^i"'+' = V(j-m)(j + l + m)^r 

where we have labeled the SU (2) index values by + and — to denote the eigenvalue. 
Combining these yield the equation 

(A + i + l)(A-i) = 

so the masses are given by M = + and M = + Explicitly, the M = + 
eigenstates are 



V2J + 1 



V'- — 7^ — y 77— — r"^' m+l 



^^In this case, we are combining the spin j representations with the spin i representation, and 
again, the eigenstates will simply be the irreducible representations. 
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where we may have < j < A'" — 1 and —j — l<m<j. The M = + ^) eigenstates 
are 



ix 



'N 



' 2j + l 



N 



Ij + l + m 



2j 



j m+l 



where 1 < i < iV — 1 and —j<iTL<j — l. 

In terms of the oscillators r] and x, the quadratic Lagrangian for the fermions 
becomes 



N-i 



t ^ t 



Here, rj and x carry fundamental and antifundamental SU{4) indices respectively, 
which we have suppressed. 

Summary of oscillator masses for the irreducible vacuum 

The spectrum of oscillators about the single membrane vacuum is summarized in table 
1, with the masses given in units of /i. 



Type 


Label 


Mass 


Spins 


50(6) X 50(3) Rep 


Degeneracy 


50(6) 


jm 




< j < - 1 


(6,2j + l) 


6(2j + 1) 


50(3) 


Oijin 
Pjm 


3^3 
J 
3 


< J < - 2 
1< j < 


(l,2j + l) 
(l,2i + l) 


(2j + 1) 
(2j + 1) 


Fermions 


Xjm 
Vl jm 


4^3 
12 ~ 3 


k<J<N-l 


(4,2i + l) 
(4,2j + l) 


4(2j + 1) 
4(2j + l) 



Table 1: Oscillator masses for the irreducible vacuum 



It is easy to check that the sum of the boson masses and the sum of the fermion 
masses are both |A"(8A"^ + 1); ^ required for the vanishing of the vacuum state zero 
point energy. 

A nice feature of the spectrum we have calculated is that the oscillator masses 
are completely independent of A^. As one might hope from the discretized membrane 
picture, N acts exactly as a short distance cutoff, providing an upper bound on the 
allowed values of the spin j (related to the wavelength of modes on the sphere). In 
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particular, the spectrum has a well defined limit for large N, which should reproduce 
the spectrum of the quadratic part of the continuum supermembrane field theory for 
this background. 

5.3 Oscillator masses in a general reducible representation 

A general classical vacuum state of the theory corresponds to an N dimensional re- 
ducible representation of SU{2), described by a block-diagonal matrix 



( Jl \ 



\ JkI 



where J/ are generators of the A^pdimensional irreducible representation of SU (2) and 
+ ... + Nk = N. 

To deal with these states, it is convenient to break up our matrices into blocks 
corresponding to the individual irreducible representations, so for example 

^kl = ^klJi 

where 1 < k,l < K are indices labeling the blocks. Then the quadratic actions spht 
up into a sum of terms which are of same form as for the irreducible representation 
above, 

= \{\f^^{{xti)\-xti + ro{XoX'^,))) 

but we now define 

Thus, whereas the matrices in the irreducible case were in the tensor product of two 
N dimensional representations of SU{2) and spht up into the representations 1 © 3 ® 
■ ■ ■ © 2N — 1, a general rectangular block Bki lies in the tensor product of the Nk and 
Ni dimensional representation and may be decomposed into irreducible representations 
with spins ^\Nk — Ni\ < j < |(iVfc -l- Ni) — 1. Thus, for example, we may write 

^kl ~ Z^K-^kUjm-i^ jm 

where F's are x Ni matrices which are eigenstates of J and J3 acting as in (5.3). 
Using this decomposition, the diagonalization proceeds exactly as above. 

^^It is important to note that for any finite /j, there will be some energy above which the perturbation 
theory breaks down and the spectrum can no longer be trusted. We discuss this further below. 
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Type 


Label 


Mass 


Spins 


50(6) X 50(3) 


,50(6) 






iWk 


-Ni 


<j<UNk + Ni)-l 




(6,2j + l) 


,50(3) 




3 ~ 3 




Ni\- 


l<j<UNk + Ni)- 


2 


(l,2i + l) 




PIT 


i 

3 




-Ni\ 


+ l<j<UNk + Ni) 




(l,2i + l) 


Fermions 


1 jm 

Xkl 


4^3 




Ni\- 


i < J < |(iVfc + iv,) - 


3 
2 


(4,2j + l) 




jm 

Vi ki 


12 ~ 3 




Ni\ + 


'^<J<l{Nk + Ni)- 


1 
2 


(4,2j + l) 



Table 2: Oscillator masses for reducible vacua 



The spectrum of oscillators about a general vacuum are given in table 2, again with 
masses in units of As before, we will have zero mass oscillators from the matrices 
corresponding to each of the non-trivial gauge orbit directions, and these will take the 
form 

Y' = [J\A] 

for arbitrary A. The number of nontrivial gauge orbit directions will be N"^ minus the 
dimension of the space of matrices that commute with J* , or 

where Mj is the multiplicity of the i'th irreducible representation. 
5.4 Physical states 

We have now computed the spectrum of oscillators about an arbitrary vacuum of the 
theory. In the fi = oo limit, one can view these classical solutions as super selection 
sectors, and the physical spectrum in each sector is obtained by acting on a Fock space 
vacuum with arbitrary numbers of creation operators for the oscillators found above. 
The only complication is that we must ensure that the Gauss law constraints, arising 
from the equation of motion for ^40, are satisfied. These A^^ — 1 conditions read 

t[X\ W] +t[X% W] + 2iIji^iIj^^" = 

The corresponding operator is ordered so that it becomes the generator of SU{N) 
transformations, and the quantum constraint requires that this operator annihilates 
physical states. In other words, physical states must be SU (N) invariant. 
Expanding the constraint about X* = ^ J\ we find 

(^i[J\ W] + z[Y\ W] + n«] + 2V^,,V't/a^ = 

In the limit — > oo, it is useful to split these — 1 constraints into two groups de- 
pending on whether or not the first term vanishes. Recalling that normalized oscillators 
in the 50(3) directions satisfy an eigenvalue equation 
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with = A^/3, and commuting this equation with J*, we find 

X[J\ m] = . 

Hence, the first term in the Gauss law constraint is non-zero precisely for the momenta 
associated with zero-mode (A = 0) oscillators = [J*, A] corresponding to gauge 
orbit directions. For these A^^ — Mf constraints (where Mi is the multiplicity of the 
ith irreducible SU{2) representation), we can ignore all but the first term in the large 
II limit, and the condition is simply that the wavefunction does not depend on the 
coordinates corresponding to gauge orbit directions (i.e. that the zero-modes are non- 
physical). 

The remaining number of constraints, J2 — 1 is exactly the number of generators 
in the subgroup of SU (N) under which the vacuum state is invariant. These constraints 
demand that physical states are invariant under this subgroup. In a vacuum with Mj 
copies of the iVj dimensional irreducible representation, the unbroken symmetry group 
will be [/(Ml) X ■ • ■ X f/(M„) (including the diagonal U(l) whose Gauss law constraint 
is trivial). From table 2, we see that the oscillator spectrum corresponding to a given 
block depends only on the dimension of the two representations defining the block, so 
for each pair of representations Aj, Nj we will have MiMj identical copies of the same 
spectrum, which form a matrix in the bifundamental representation of U{Mj) x U{Mj) 
(or in the adjoint of U{Mi) for the case Ni = Nj). The Gauss law constraint is then 
satisfied simply by forming U{Mi) x ■ ■ ■ x U{Mn) invariants out of the matrix creation 
operators. Thus, a general physical state will be obtained by acting with an arbitrary 
number of traces of products of these matrix creation operators on the fock space 
vacuum. 

We have now specified the complete spectrum of physical states of the matrix model 
in the // — > oo limit. 

5.5 BPS States 

We have argued in section 4.1 that the unusual supersymmetry algebra of the ma- 
trix model allows the possibility of states carrying angular momentum which preserve 
various fractions of the 16 supersymmetries generated by Q. Now that we have the 
complete spectrum of the theory in the /i — oo limit, it will be possible to explicitly 
identify all BPS states. 

To begin, it is convenient to write the anticommutation relation for the Qs in the 
original SO {9) notation, 

{g«, Qp} = 25^f,H - ^(^'''YXpM'' + ^(l'''r%pM'^' (5.9) 

We may choose to label states by their eigenvalues for a set of Cartan generators of 
the S0{3) X S0{6) algebra, which we take to be M^^, M^^ and M^^. The 7 

matrices coupling to each of these generators all commute with one another and may 
be represented by 

73 ^ (7^x1x1x1 7^23^^ = 1x0-^x1x1 
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^12367 ^ ixlxa^xl 7^23^^ = -1 X X (T^ X 1 . (5.10) 



where the last definition is imphed by 

„,12345^,12367^,12389 ^,123456789 i 

7 7 7 = —7 = — i . 

Any non-trivial BPS state must have non-zero angular momentum and therefore will 
lie in some non-trivial representation of 5*0(3) x S0{6). For this 5*0(3) x 5*0(6) 
multiplet, there will be a unique highest weight state and for this state, we will have 
(M^-^) = for all of the non-Cartan generators of the algebra. Using the representation 
(5.10), the expectation value of (5.9) on this state is then 

{^lj\{Qa:Qf3}\^)=2SapAp 

where the sixteen diagonal elements are given by two copies of the values 

A = H + 6i^Mi2 + e2^M'' + esf^M^' - e^e^^M^' 
3 6 6 6 

and €i — ±1 may be chosen independently. Highest weight BPS states will then be 
exactly those for which some of these eigenvalues vanish. It is now straightforward to 
describe the complete set of BPS states in the spectrum by inspecting table 2. 

The eigenvalues of 12{H/ fi, M^^S, M^^G, M^'^ /6, M^^ /6) for the individual oscil- 
lators (suppressing the block indices) are 

Pj^: (4i, 4m, 0,0,0) 
0(jm- (4j +4,4m,0,0,0) 

Xjm ■ (4j + 2, 4m, ±2, 0, 0), (4j + 2, 4m, 0, ±2, 0), (4j + 2, 4m, 0, 0, ±2) 
Xjm-. (4j + 3,4m,l,l,l), (4j + 3, 4m, 1,-1,-1) 

(4j + 3,4m, -1, 1, -1), (4j + 3,4m, -1, -1, 1) 
rjjm- (4j + 1,4m, -1,-1,-1), (4j + 1, 4m, 1, 1, -1) 

(4j + 1, 4m, 1, -1, 1), (4j + 1, 4m, -1, 1, 1) 

The {Q, Q} eigenvalues are obtained by dotting these vectors with 

— (^1, ei, 62, 63, —6263; , 

and a state corresponding to one of these oscillators is invariant under the pair of super- 
charges corresponding to (ei, 62, 63) whenever the result is zero. It is straightforward to 
check that this dot product is always greater than or equal to zero (guaranteed by the 
BPS bound), and that it may only be zero for m = ±j. For this choice, we find that 
each P, T], X, x, and a oscillator commutes with 4,3,2,1, and pairs of supercharges, 
respectively. 

^•^Note that the mode a;oo is an exceptional case and the number of pairs of supercharges preserved 
is 4. 
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The oscillators which preserve each supercharge are given by 



+ + + : 


(^jj ' 


^4-i5 
•^jj ' 


^6-i7 
•^jj ' 


^8+j9 
•^jj ' 










+ + - : 


(^jj ' 


4-j5 
•^Jj ' 


•^Jj ' 


8-i9 










H h : 


(^jj ' 


•^ii ' 


6-j7 

■^ii ' 












H : 






^6+i7 




+++ 









Here, we have labeled fermion oscillators by the signs of their (M"^^, M^^, M*^) eigen- 
values and denoted + ix^ by x^"*"'^. The values on the left denote the sign of ei, 
62, and €3 for a given pair of supercharges. The oscillators that commute with super- 
charges corresponding to ei = — 1 are given by an identical table with the signs on the 
left reversed and m = —j for all oscillators instead of m = j. 

Since all of the quantum numbers are additive, a state obtained by acting with an 
arbitrary combination of oscillators in a single row of the list (5.11) will be invariant 
under the corresponding supercharge. Thus, up to 5*0 (3) x 5*0(6) rotations, all BPS 
states in the theory may be obtained by acting on the Fock space vacuum with an 
arbitrary combination of oscillators in any given row above. States for which all of the 
oscillators appear simultaneously in k rows of the table above will be left invariant by 2k 
supercharges. For example, a general class of 1/4 BPS states preserving 8 supercharges 
is 

n(^5.-rio)- 

where b is the creation operator associated with (3. It is interesting to note that the 
simple 1/4 BPS states, 

(4)"|0) 

correspond classically to the rotating ellipsoidal membrane state found by Bak [24]. 

We have thus given the complete set of BPS states of the matrix model for the 
^ = 00 limit. It would be interesting to determine which of these remain BPS as we 
turn on interactions (at finite n). This issue is discussed further in the remarks of 
section 8. 

6 Perturbation theory 

In this section, we move away from the /i — 00 limit and study perturbation theory 
about the free theory described in the previous section. Specifically, we would like 
to determine the range of parameters (A^, R, n) and energies for which perturbation 
theory is valid. As a first step, we rewrite the interaction Hamiltonian in terms of 
normalized oscillators to determine the value of the "bare coupling" corresponding to 
different vacua. We then perform explicit calculations of the energy shift for the first 
non-trivial excited states about the X — vacuum and the irreducible vacuum to 
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determine the actual effective coupling controlling the perturbation expansion. Since 
these states are non-BPS, we expect their energy shifts to be generic. For fixed values of 
the parameters, we then estimate the height of the potential separating different vacua. 
This corresponds to the energy above which perturbation theory becomes inapplicable. 



6.1 Expansion in terms of normalized oscillators 

In the previous section, we have explicitly diagonalized the quadratic actions, writing 
them in terms of canonically normalized oscillators. It is useful to consider expanding 
the interaction Hamiltonian in terms of these oscillators, in order to more accurately 
determine the parameter in the action suppressing interaction terms. 

In general, the cubic interactions in the matrix Lagrangian take the form 

H^r^ Tr ^[^2,^3]) 

while the quartic interactions take the form 

//4~ (^^) Tr {[A,,A2][As,A4) 

where Ai represent any of our matrix variables. 

As an example, we consider a vacuum corresponding to K copies of the Ni dimen- 
sional irreducible representation, such that KNi = N. In this discussed in 
section 5.3, it is convenient to spht the matrices A into A'^i x A^i blocks, and we find 
from the analysis of section 5.2 that the individual blocks may be written in terms of 
normalized oscillators Aj^ as 

^'^^ ~ ,y]Y^^J"*^jrn^m^ (6-1) 

where Cjm are coefficients independent of A^. To evaluate the coefficients in the cubic 
and quartic actions, we note (see, for example [21, 22]) 

where b is zero when ji + j2 + J3 is even and otherwise 

bnrn.nrn.j,m, = (-l)^2Ari ^(2ji + l)(2i2 + l)(2j3 + l) 

/ jl J2 
y mi 777.2 

Explicit formulae for the Wigner 3j and 6j symbols appearing here may be found in 
[23] . Using these, it is straightforward to show that b has a finite limit as A" — > 00 with 
corrections of order 
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Using (6.1) and (6.2) and recalling the normalization condition (5.3) for F's, we 
may now rewrite the interaction Hamiltonian in terms of normalized oscillators. We 
find schematically 

3 

/ R \ 2 

bcccAAA 





H4 ~ bbccccAAAA 

For large the leading terms in the b and c coefficients are independent of N, so 
we conclude that the coupfing corresponding to the vacuum with K spheres of radius 
r = [iiNi/QR) is given by 

For a fixed value of A^, we see that the irreducible vacuum appears to have the weakest 
coupling, Qirr — {R/{fJ,N))2 while the X — vacuum is the most strongly coupled, 

3 

with gx=o = {R/fJ')^- 

In systems with many degrees of freedom, it is often the case that the perturbation 
expansion is controlled by some effective coupling which is different than the bare 
coupling appearing in the Lagrangian. To determine the true effective coupling in our 
case, we now turn to some explicit perturbative calculations. 

6.2 X = vacuum 

The X — vacuum is invariant under the full U {N) symmetry group of the model. 
As discussed in section 5.4, it is then convenient to introduce matrix creation and 
annihilation operators transforming in the adjoint of this symmetry group. 



X' = /|(A + 4) 

X« = V3{Aa + Ai) 



in terms of which the quadratic Hamiltonian becomes 

H2 = Tr {^AlA, + ^-AiA, + ^^j^'^^i.) . 

Physical states satisfying the Gauss law constraint will then be obtained by acting with 
traces of products of these matrix creation operators on the Fock space vacuum, for 
example 

IV (44)|o) ■ 



29 



State 


,50(6) X ,50(3) Rep. 


Energy 


Degeneracy 


|0) 


(1,1) 





1 


Tr {Aim 


(6,1) 


a 

6 


6 


Tr (^^^")|0) 


(4,2) 


M 
4 


8 


Tr (4)|0) 


(1,3) 


M 

3 


3 


IV(At)Tr (4)|0) 


(1,1) + (20,1) 


M 
3 


1 + 20 


Tr (AlAl) - ^Tr (Al)Tr (4)|0) 


(1,1) + (20,1) 


3 


1 + 20 



Table 3: Lowest energy states about the X = vacuum 



The lowest lying states are given in table 3. 

Wc can now investigate the change in energy of these states upon turning on the 
interaction Hamiltonian, using perturbation theory in ^. Note that only the last state 
in the list above can get a shift in energy, since the other states are excitations in the 
U{1) part of the theory, which is free. 

It is easy to see that the energy shift at first order in perturbation theory 

will be zero for all states, since all terms in H3 have non-zero H2 eigenvalue. Thus 
H3 acting on any state gives a combination of basis states all of whose energies are 
different than that of the original state. 

At second order in perturbation theory, the energy shift is given by the usual 
formula, 

A = A4 + A3 

= {m^m + T.^r^\{n\Hsm\' (6.3) 

Before calculating the energy shift for the excited state, it is useful to verify explic- 
itly that the vacuum state receives no correction. In appendix C.l, wc show that for 
the vacuum state, A4 = — A3 = ^{N^ — N), so the vacuum state does remain at zero 
energy to second order in perturbation theory. The vacuum wavefunction does receive 
a correction at first order in perturbation theory, given in equation (C.3). 

We now compute the vacuum shift for the lowest energy non-trivial state (i.e. in- 
volving more than U{1) oscillators), 

1^) = ^=i— (TV (44) - ^Tr (4)Tr (4))|0) 
^12(A^2 _ 1) W 

This normalized state has been chosen to be orthogonal to the state 

Tr {Ai)Tr (4))|0) 
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which cannot receive any energy correction. 

The calculation of the energy shift to second order in perturbation theory is per- 
formed in appendix C.l, and we find 

Thus, the effective coupling controlling the perturbation expansion appears to be (ig- 
noring numerical factors) 

N{R/i^f = g^N 

3 

where g = {R/ jj)^ was the bare coupling suppressing the interaction terms in the action 
near X = Q. 

6.3 Perturbation theory about irreducible vacuum 

Wc next consider perturbation theory for the single membrane vacuum, corresponding 
to the irreducible representation of SU (2). Here, the low-energy states are obtained by 
acting with creation operators corresponding to the oscillators in table 1. The states 
with lowest energy are summarized in table 4, where we have defined a^^, a^-^, and 
hjm to be canonical annihilation operators associated with the oscillators x"^, ajm, and 
(5jm respectively. 



State 


SO{i6) 


X 50(3) Rep. 


Energy 


Degeneracy 


|0) 




(1,1) 





1 


Ko)^|o) 




(6,1) 


M 
6 


6 






(4,2) 


4 


8 


&Lio) 




(1,3) 


M 

3 


3 


Ko)^(«oo)^|0) 


(1, 


1) + (20, 1) 


M 
3 


1 + 20 


4o|0) 




(1,1) 


3 


1 



Table 4: Lowest energy states about the irreducible vacuum 



Note that Xqq, f3im and rji^ correspond to oscillators in the free U{1) part of the 
theory, so it is only the last state that can receive an energy shift in perturbation 
theory. This is the breathing mode of the sphere. 

In appendix C, we compute the leading energy shift for an arbitrary occupation 
number of this breathing mode, i.e., for the state 

4t(«oo)"|0) . 
\/n\ 
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The result, coming at second order in perturbation theory from the terms (6.3) is 



(Note that for n = 1 the energy shift A is zero.) In this case, we see that the leading 
correction for low energy states is controlled by the "bare" coupling without any 
additional A^-dependence and is therefore finite in the large N limit. 

6.4 Validity of perturbation theory 

In the previous sections, we have computed the leading correction to the energies of 
the first nontrivial excited states above the irreducible single-membrane vacuum and 
the X = vacuum. We found that the effective coupling controlling the perturbative 
corrections is given by A = g^K^ where K — \ and gf^ = (^)2 for the single membrane 

case and K — N and g"^ — (^)i for the X — Q vacuum. Note that K corresponds 
to the number of coincident membranes characterizing the vacuum and in this sense A 
behaves like a 't Hooft coupling. 

In situations where perturbative corrections are small for low lying states, we would 
now like to estimate the energy above which perturbation theory breaks down. In terms 
of the classical potential, we expect this to occur when we move far enough away from 
the minimum so that the interaction terms in the potential become of the same order 
of magnitude as the quadratic terms. 

To obtain an initial, crude estimate of this energy scale, consider starting at this 
irreducible vacuum and exciting the breathing mode considered in the previous section, 
that is, take y = J*. Plugging this in to the classical potential, we find 



Thus, we see that the cubic and quartic terms in the potential will be small compared 
to the quadratic term when a -C 1. In terms of the energy, this translates to the 
condition (ignoring numerical factors) 



Our crude estimate suggests that the quadratic approximation to the potential is good 
for energies less than ^x/g^. 

^**The other vacua are characterized by a couplings intermediate between these values. 




V 
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Figure 2: Low energy path between nearby vacua 



We expect this estimate to be true for a generic direction away from the minimum 
However, it turns out that for certain paths away from the minimum, the interaction 
terms can become important at a much lower energy. In appendix D, we provide an 
explicit path between the irreducible vacuum labeled by N and the nearby vacuum 
labeled by (A^ — 1, 1) along which the maximum energy is of order {/jlN/RYjj: — 
rather than l/g"^}^ This path corresponds to the discretized version of a continuum 
picture in which a thin tube of membrane extends from the south pole of the sphere to 
the origin where a small second sphere can grow from the tip of the tube, as depicted 
in figure 2}^ 

Thus, our crude estimate for the energy at which interaction terms become impor- 
tant was in fact too large by a factor of N. At least along a specific direction away 
from the minimum, we find that interactions must be important at an energy 



From this result, one might conclude that the condition for a perturbative regime to 
exist should be obtained by demanding that this energy is significantly greater than 
the energy of the lowest lying modes, that is, g'^N <S 1. On the other hand, we 
have seen explicitly that the leading perturbative correction for certain low energy 
states was small as long as (7^ <^ 1, a much weaker condition. It is possible that the 
low energy paths between vacua we have found are very narrow for large (a set of 
very small measure in the space of paths away from the minimum) and therefore that 
perturbation theory may be valid for larger energies than the height of these paths 
would indicate. In particular, from the calculation of appendix C.2, we see that the 
the leading perturbative correction to the breathing mode with occupation number n 
is small when 71^ g'^ <^ 1, in other words when E -C iJ./g. Further work will be required 
to decisively determine which of these conditions (or some intermediate condition) 
actually governs the validity of perturbation theory. 

^^In terms of the natural metric Tr ((X"^)^) on the space of matrices, this choice for the final state 

represents the closest vacuum. 

^^We thank Washington Taylor for a valuable discussion about this. 
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7 M-theory limit 



In this section, we discuss the imphcations of our results for the M-theory hmit of the 
matrix model. 

7.1 Implications of perturbative results 

M-thcory on the pp-wave background should be described using the matrix model by 
taking the .'V ^ cxo limit with fixed = N/R and fixed fi. The resulting theory 
is characterized by the single boost- invariant dimcnsionlcss parameter np'^l^. In this 
background, a single membrane with momentum pi in the pp-wave background will 
expand into a spherical BPS "giant graviton" configuration with radius ixpf [6]. In 
general, the total momentum p'^ may be divided among any number of membranes, 
and each of these membranes will have a radius proportional to its momentum. This 
is precisely the picture we get from the large N limit of the matrix model. 

In section 5, we calculated the spectrum of the matrix model in the limit of large 
ji. We found that the individual oscillator masses were independent of N. For the 
irreducible vacuum, acts as a cutoff controlling the maximum allowed frequency of 
the modes. Thus, in the large limit, we find a well- defined spectrum, which we may 
associate with the spectrum of the free part of the field theory on a single spherical 
supermembrane in the pp-wave background. 

For reducible representations involving k copies of a single irreducible representa- 
tion, we obtain k"^ copies of the spectrum for the corresponding irreducible represen- 
tation, suggesting a non-abelian structure for coincident membranes. For reducible 
vacua involving different irreducible representations of dimension A'^i and iV2, corre- 
sponding to non-coincident membranes, we found that the spectrum corresponding to 
off-diagonal degrees of freedom connecting the two representations has a lower bound 
on oscillator masses proportional to ||iVi — iV2|. Thus, if we fix the radius parameters 
(fiNi/R) of the two membranes and take the large limit, we find that all of the 
off-diagonal degrees of freedom develop infinite masses, going like N{pf —p'D/p'^- 

It is not clear whether the M-theory limit is compatible with a valid perturbation 
expansion. In explicit perturbative calculations, we found that the leading corrections 
to the energies of low-lying states above the singlc-mcmbranc vacuum are finite in the 
the M-theory limit and small when ^p'^ ^ 1. On the other hand, we have shown in 
appendix D that there exist paths in configuration space between any two vacua with 
maximum energy of order n^np'^Y /N , which goes to zero in the M-theory limit. This 
indicates that in the large N limit, the classical vacua are no longer isolated, but rather 
are all connected by finite, zero-energy fiat directions in the potential. Physically, these 
fiat directions correspond to decreasing the radius (and light-cone momentum) of any 
membrane while increasing the radius of any other membrane at the same rate (or 
growing a new membrane at the origin). Quantum mechanically, these fiat directions 
would seem to spoil any perturbative study of the model, since none of the classical 
vacua are long-lived and the low-energy wave-functions would presumably spread out 
over all the classical vacua. 
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7.2 The supermembrane field theory 

From the supermembrane point of view, the matrix model can be thought of as a 
regulated version of the light-cone supermembrane field theory on some fixed topology 
worldvolume. When interpreted in this context, our results have a number of interesting 
implications. 

To be definite, consider starting from the supermembrane Hamiltonian (2.29), and 
choose a spherical worldvolume of radius 1 labeled by coordinates Xi,X2,X3 with xf + 
X2 + x^ = l}"^ The potential term in the bosonic part of the Hamiltonian is 

V oc {^''''{X\X^} - ^y^^^j + {X\X^Y + 2{X\X''Y + (^)'(X")' . (7.1) 

where {A, 5} = e^'^^XidjAdkB is the canonical Poisson bracket. This potential has a 
minimum when X" = and 

{X^X^■} = e^^■*^^X*^ (7.2) 
3 

In the continuum case, there are only two vacua corresponding to non-singular field 
configurations, namely X — and X^ — Rsph^^ where Rsph = HP'^ /Q- 

By analogy with the matrix model analysis above, one might have been interested to 
study the field theory expanded about the large membrane solution. The resulting the- 
ory has an unusual gauge symmetry coming from the symmetry under area-preserving 
diffeomorphisms , 

5Y' = -Rsphe^'^x^d^K + {Y\ A} 

One may attempt to decompose Y into fields with more standard gauge transformation 
properties as 

= --^-x-y, A^-^xxY, (7.3) 

^sph ^sph 

To leading order, the fields A and have the usual transformation properties of a gauge 
field and a scalar, and furthermore A appears in the action as a massless gauge field 
while (f) is massive. However, there are also nonlinear terms in the gauge transformation 
which may only be removed by adding higher order terms in the definition of A and 
0. These introduce even higher order terms in the gauge transformation, and so forth, 
indicating that the definition of a proper gauge field with transformation property 
5Ai = diA — {Ai. A} would require an infinite scries of terms in Y. Thus, for now we 
shall consider the theory written in terms of Y. 

At first sight, this theory looks non-renormalizable, since the { , }^ terms have 
dimension 6, suggesting for example that one would find divergent one-loop contribu- 
tions in the quadratic effective action going like A^. In the matrix regularized theory, 
the role of the UV cutoff is played by A^, so we might have expected A^^ behavior 
in the mass shift we have calculated for the breathing mode. In reality, our result is 



^^In terms of the usual spherical coordinates (j),9, we have xs = cos 6, xi = sin cos 0, and X2 
sin ^ sin ^ and the Poisson bracket on the sphere is defined by {6, (j)} = 
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0{1) in the large N limit, indicating that the theory may have much better renormal- 
izability properties than it appears to. This is at least partly due to supersymmetric 
cancellations. Another possible explanation is that we should really consider the the- 
ory as a ^ — > limit of a noncommutative 2+1 dimensional field theory on a sphere 
with only quartic interactions. Such a theory would be expected to give only linear 
divergences in the quadratic efi^ective action, or even weaker behaviour in the super- 
symmetric case. It would be quite interesting if the Poisson bracket field theory turned 
out to be perturbatively renormalizable. 

We noted that the supermembrane field theory has two nonsingular classical vac- 
uum solutions, however, from the matrix theory point of view, we have seen that the 
theory has a continuum of other vacua connected to these by finite zero-energy fiat 
directions. In the continuum theory, these vacua will appear as singular field configu- 
rations. To understand what these look hke, recall that one path between the single 
membrane vacuum and a vacuum with two spheres involved forming a thin tube of 
membrane from the south pole of the original sphere which moves up to the origin and 
nucleates another sphere. From the membrane field theory point of view, this corre- 
sponds to nucleating a small bubble of domain wall at some point on the sphere, where 
the inside of the domain wall maps to the new sphere, and the domain wall itself maps 
to the tube. In the final vacuum state with two spheres, the tube has shrunk to zero 
size, so the domain wall becomes singular. Thus, the fiat directions we have found cor- 
respond to the existence of massless singular domain walls, or tensionless strings, in the 
supermembrane field theory. Presumably, these nonperturbative massless excitations 
make the theory extremely difficult to study. 

8 Summary and Discussion 

In this paper, we have studied the matrix model proposed by Berenstein, Maldacena, 
and Nastase to describe M-thcory on the maximally supersymmetric pp wave solution 
of eleven-dimensional supergravity. We have seen that the model may be derived di- 
rectly as a regularized form of the light-cone supermembrane action in the pp-wave 
background, or alternatively (with certain assumptions) from DO-branes in a corre- 
sponding type IIA background. 

We noted that for large values of the mass parameter fi, the theory may be ex- 
panded about its classical vacua to give a quadratic theory with interactions suppressed 
by powers of If/i. In the large /j, limit for fixed N and R, we determined the complete 
physical spectrum of the theory and found that the individual energy levels are rational 
multiples of fi independent of A^. We characterized the complete set of BPS states in 
this limit, finding infinite towers of states preserving 2,4,6, and 8 of the 32 supercharges 
(in addition to the 1/2 BPS vacuum states). We then considered perturbation theory 
about the vacua, and found that vacua corresponding to more spheres at smaller radii 
are generally more strongly coupled. For the most weakly coupled vacuum, correspond- 
ing to a single membrane sphere, we found that the leading corrections to the energies 
of low- lying excited states were small when g'^ — l/(//p+)^ -C 1, independent of N. 
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On the other hand, we saw that the various vacua are connected by paths along which 
the maximum energy is {np'^Y /N indicating a more restricted regime of vahdity for 
perturbation theory incompatible with the M-theory limit. 

A very interesting direction to pursue would be to establish which of the BPS states 
appearing in the fi = oo limit remain supersymmetric for finite fi. Typically, BPS mul- 
tiplets may only become non-BPS by combining with other BPS multiplets to form 
larger non-super symmetric representations. By better understanding the representa- 
tion theory of the pp-wave superalgebra, it may be possible to show that certain BPS 
multiplets arc protected. The existence of protected states would provide reliable in- 
formation about the spectrum even in the regime where perturbation theory is invalid. 
Another approach to determining which states remain BPS would be to perform a 
direct calculation of the energy shift in perturbation theory, as we have done in this 
paper for two simple non-BPS states. Since the energy of BPS states is determined 
exactly by their angular momenta, any state which remains BPS when interactions arc 
turned on must have zero energy shift in perturbation theory. We note that at least 
one of the 1/4 BPS states we have found at = oo should remain BPS, namely the 
rotating ellipsoidal membrane solution considered in [24]. 

We would now like to briefly address the issue of vacuum stability. We have seen 
that classically, and to leading orders in perturbation theory, the various vacua are all 
degenerate and have zero energy. It is not clear whether this will be true nonpertur- 
batively in the full quantum theory. We know that for = there exists precisely 
one normalizable zero-energy state, however it is possible that other zero-energy states 
(corresponding to multiple gravitons) exist which become non-normalizable as // — > 0. 
If the classical vacua do not remain degenerate, then the exact wavefunctions of the 
theory will presumably not be localized around any one of the classical vacua, and tran- 
sitions will occur. On the other hand, if fi is large, the low-energy wavefunctions we 
have been discussing should serve as useful approximate wavefunctions over time scales 
much shorter than the typical decay rate from one vacuum to another. We expect that 
for large values of the path integral governing the transition rate between different 
minima of the classical potential (corresponding to membrane transition amplitudes) 
will be dominated by classical instanton solutions so that the rates may be estimated 
by a stationary phase approximation. Some examples of these instanton solutions have 
been worked out in [25]. 

It is interesting to note that the more commonly studied case of type IIB string 
theory on the maximally supersymmetric pp-wave solution of type IIB supergravity 
should behave very much like the M-thcory case when the string coupling is taken to 
be of order 1. In particular, along the range of f^p^a' with fixed gg = 1, the classical 
vacua of the theory are characterized by collections of D3-brane giant gravitons with 
the sum of the squared radii proportional to np'^ (in units where a' = 1). It would 
be interesting to determine whether a matrix model for the DLCQ of type IIB string 
theory on this background (discussed recently in [26]) also contains a perturbative 
regime for large jj,. 
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A Fermions in SU{2) x SU{4) notation. 

The massive matrix model preserves an 5*0(3) x 5*0(6) ~ SU{2) x 5'[/(4) subgroup of 
the original 50(9) symmetry of the BFSS model. It is convenient to write the action in 
a form for which the 5'0(3) x 5'0(6) transformation properties of all fields are manifest. 
This is already the case for the bosons, so we turn to the fermionic fields. Under the 
decomposition 

50(9) ^ 50(6) X 50(3) 
the spinor representation splits up as 

16 ^ (4, 2) + (4, 2) 

where we take / as a fundamental SU (4) index and a as a fundamental SU (2) index. 
Our spinors obey a reality condition, which in the reduced notation becomes simply 

To decompose the terms in the action, we introduce the matrices gjj which relate the 
antisymmetric product of two 4 representations of SU (4) to the vector of 50(6). Then, 
we may choose"*^^ 

and 

/ lxg"\ y^f-^^xl \ 

' \^ 1 X (g«)t ) ' \^ a'xl ) 

These Dirac matrices satisfy the Chfford algebra as long as we take normalizations so 
that 

^^In this representation, the Dirac matrices and fermions arc not real, however the spinors still obey 
a reality condition which reduces them to 16 independent real components. 
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We may then rewrite the quadratic fermion terms in our original action as^^ 
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B Symmetry algebra in terms of matrix generators. 

In this appendix, wc review the symmetry algebra of the matrix model and provide 
explicit expressions for the generators in terms of the matrix model variables. 

The bosonic generators include the Hamiltonian H, the light-cone momentum P+ 
(a central term in the algebra) and rotation generators of the SO (3) x SO (6) symmetry 
which we denote by M^^ and M"-^ respectively. In addition, there are generators Oa and 
Oj with nonvanishing commutation relations 



Q^o) Oft 



= P+S. 



ab 



6 

iSacO-b - iSbctta 



These are the creation and annihilation operators corresponding to the decoupled U{1) 
part of the theory which describes a particle (the center of mass coordinate) in a 
harmonic potential. 

The algebra has in addition 32 fermionic generators. The first 16, labeled by qia 
change the overall polarization state in the C/(l) part of the theory, while the other 16, 
labeled by Qia act non-trivially. 

These generators obey anticommutation relations 



{Q^''',Qjfs} 

{Qia, Qjp} 
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Since we are working with fermions that are not exphcitly real, we need to replace with ^FT. 
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2 



{q'"'',Qjp} = - 

and have nonvanishing commutation relations with the bosonic generators given by 



[H, qia 

[H, Qia 

Qiai 



-qio 



Q 



la 



ii 

12 



-^^|lg'iJ^a,q''' 




^CTa qip 



with additional commutation relations between the fermionic generators and M's ap- 
propriate to the 5*0(3) X 5*0(6) transformation properties of the qs and Qs. 
These generators are realized explicitly in the matrix model as 



qia 

Qia 



R 



IV (1) 



1 

Vr 



Tr 



12R 



I3R 
2/1 




in") 



= Tr (X^n^' - X^n^ + ie'^^il)'^ a^^lj) 
= Tr (X"tf - X^n" + -t/'^g^V) 



\r' 



RTy ( {W - z^X«)g?^e«^V'tJ/3 _ (n^ + , Ax>;/3^,^ 



ab\J„ 



where g"^ = |(g"g^^ ~ g^'g^")- The Hamiltonian H may be read off from the action in 
section 4.1. Using the quantum commutation relations 



1 1 

mn\ 



iS''''SknSlm {V'ir> (^J/3)mn} = ^J^'^hn^lm 



]Ia 



I xa} 



we may verify explicitly that the algebra given above is satisfied. 

For convenience, we also provide the supersymmetry algebra in the original 5*0(9) 
notation of BMN, 

{Qa,Qp} = 25.^//-|(y^V)a/3M^^ + |(y2V)a/jM«'' 
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where 



C Pert ur bat ive calculations 

In this appendix, we perform exphcitly the calculation of the energy shift for the lowest 
non-trivial excited states above the X — Q vacuum and the irreducible vacuum. 



C.l X = vacuum 

Writing out the cubic interaction in terms of the matrix creation and annihilation 
operators, we find 

3 / I — 

= f-Vxr \\ he'^%A,A^Ak + ZA\A^Ak + ZA\A]Ak + A\A]Ai) 



For the quartic term, it will be more convenient to keep the coordinate notation 



It is easy to see that the energy shift at first order in perturbation theory 

will be zero for all states, since all terms in have non-zero H2 eigenvalue. Thus 
acting on any state gives a combination of basis state all of whose energies are different 
than that of the original state. 
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At second order in perturbation theory, the energy shift is given by the usual 
formula, 

A = A4 + A3 

= (V'IW) + ETr-^IHW)r (c.i) 

Wc start by verifying that the vacuum state stays at zero energy. The contribution 
to the vacuum shift from is given by 

A4 = {0\H^\0) 
27 

= — (A^^-A^) + 135(A^^-A^) + 81(A^^-A^) 

where the three terms come from the ij, ab, and ai terms in the quartic potential 
respectively. 

To compute A3, note that 

/3 
-? 
8 



i/alO) = ^-i^l8iN^-N)\A) + V3^A8{N^ - N)\B) (C.2) 
where we have defined orthonormal states 



\B) = 48(jV3_jv) '^<-'g^1^(^">'-"'-^i)l°) 

Thus, the contribution to the vacuum shift from H3 is 

A3 = -^mHs\0)\' - ^\{B\Hsm' 
= -'^{N^ - N) - 216{N^ - N) 



-A. 



4 



Thus, we have shown that the vacuum state remains at zero energy to second order 
in perturbation theory. Using (C.2), we may also compute the correction to the vacuum 
wavefunction to first order in perturbation theory, and we find the result 



\vac) = (10 > -J-iXJl8(N^-N)\A) - 18XJ(m - N)\B)) 

' ' ^i + Mip(^3_Ar)^' V8 V ^ ^1 / J\ n 

We now compute the vacuum shift for the lowest energy non-trivial state (i.e. in- 
volving more than U{1) oscillators), 
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This normalized state has been chosen to be orthogonal to the state 



TV (Al)Tr (Aim 

which cannot receive any energy correction. 

To compute A3, we note that may be written as a sum of terms with specific 
H2 eigenvalues, 

such that 

where H2 is the quadratic Hamiltonian and all the E^s are distinct. 
Using this decomposition, we find 

A3 = j:^r^\{n\Hsm' 

= E-;^(^l(^3)V)Hi^3l^) 

n,i,j * 

= j:-U^\iHmn){nm^) 

n,i * 

Thus, we may write A as 

A = (v^|i^4|V') -E;^(V'I(^3)^^^IV^) (C.4) 

i ' 

Since the vacuum energy shift was zero, we can ignore all "disconnected" con- 
tributions in which the creation operators from the initial state contract only with 
annihilation operators from the final state since these terms will be the same as for the 
vacuum shift. A further simplification arises from the fact that the interaction vertices 
are written in terms of commutators of X's with the result that Tr (A' ) or Tr (A) from 
the initial or final state contracted with any of the interaction terms will vanish. Thus, 
we may ignore the 

-^Tr (4)Tr (4) 

piece since it contributes only to the disconnected contributions. 
The connected contributions to A3 come from the terms 



Hi = V3e^f,g'}jTr{^^'^'^^^'Al) 
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with H2 eigenvalues £^1 = |, and E2 = ^. The two contributions to A3 are then 



-1447V 



and 



= -72 AT 

The connected part of A4 also receives two contributions, coming from 
Hi = -^Tr {X'^X^X^'X^ - X'^X^'X^X^) 

and 

Hi = -Tr {X^X^X^X' - X^X^X'X') 

The computation may be carried out efficiently using diagrams with double-line nota- 
tion, making use of the propagators 



Evaluating the various figure eight diagrams we find 

Al = - i2(iV^2 _ 1) ■ ^(T^ (-4c^c)TV iX'^X'X'^X' - X'^X'^X'X')Tr (44)) 
= 2707V 



and 



A^ = ~ i2(jV2 _ 1) ^^^ i^cAc)Tr {X'^X'X'^X' - X'^X'^X'X')Tr (44)) 
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Thus, we find that the total energy shift for the lowest nontrivial excited state IV') is 
given by 



3 

= 108iV/i l^^j 

Thus, the effective coupling controlling the perturbation expansion appears to be (ig- 
noring numerical factors) 

Xx=o = N{R/fj,f = g^N 

where g = {R/ jj)^ was the bare coupling suppressing the interaction terms in the action 
near X — Q. 



C.2 Irreducible vacuum 

We now compute the energy shift for an arbitrary occupation number of the breathing 
mode of the membrane sphere for the irreducible vacuum, 

\^) = 4^(4o)lO) . 
\/n\ 

We want to calculate the energy shift upon turning on the interaction terms 

= (-)iTr (--[J\ X"] - -[J\ Y^\\Y\ + -e^^^Y'Y^Y^ 

^ \ 3 3 3 

H, = (^fTr (-l[y^y^f -i[x",x'']2-l[n^1') (C.5) 

As in our calculation for the irreducible vacuum, the leading correction to the energy 
will come at second order in perturbation theory and will be given by 

I^^ {^H,\ij) -Y^^mHiyHm (C.6) 

i * 

where, as discussed in the previous section, we have expanded 
such that 

where H2 is the quadratic Hamiltonian and all the EiS are distinct. 

Again, we need only evaluate the connected contributions, since the disconnected 
pieces will be the same as for the vacuum state and therefore cancel. For all connected 
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contributions, the as and a^s from the initial and final state must contract with some 
Ys in i?3 or H4, and these Wick contractions may be performed using the relations 



i'^^mn) ~ O^mn'^^) ~ \l at t at2 T\'^mn 



N{N^ - 1 

In other words, Wick contracting a Y with a or replaces it with a J. 
The connected contributions may be divided as 

A = A + A 

where A denotes the contributions for which all but one of the as from the final state 
contract with a^s from the initial state, while A denotes those for which all but two of 
the as from the final state contract with a^s from the initial state. Taking into account 
the normalization factor for the state and the combinatorical factors associated with 
pairing as from the final state with a^s from the initial state, we find 

A = n • A{n = 1} 

A = (^l^A{n^2} (C.7) 

We begin by calculating A{n = 1}. 

The term in A will involve terms in the quartic action with two F's replaced 
by J's, while the H^H^ term in A will involve terms in the cubic action with one Y 
replaced by a J. The calculation therefore reduces to evaluating correlation functions 
involving the following operators which are quadratic in the fields. 

^ = Tr {ie'^^f[Y^,Y^]) 
B = Tr ([J^y^f) 
C = Tr([/,X"]2) 
D = Tr (V'V^[J^V']+V'(^iJ^V'1) 
In terms of the normalized coordinates, these are given by 

c = -Ej(j + i)(4J^4m 

t t t t 

^ ~ ~ 2^iVjmVjm ~ VjmVjm) "I" "I" 2^iXjmXjm ~ XjmXjm) 

The expectation values of these operators may be calculated to be 

3(iV-l) 



(A) 



2N 



(B) = -^{N - 1){4:N^ + AN - 3) 

(C) = -6N{N^-1) 
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using the propagators 



t 

{{Oijm) Oi~jrh) = ^33^'^^'^2{i + 1) 

3^ 
27 



The contribution may now be quickly evaluated. We find 

A^^^ = (a(-lTr ([XM^f ))at) = • 2 • ^^^^ _ (C) 

= 36 

while 

= 12-9- 



N'^{N + l) 

Note that in the first step we have used the fact that 

Tr {[J\ Y^\[Y\ J^]) = Tr ([J^ Y'f) + Tr {[J\ J^\\Y\ Y^\) 

but the first term on the right involves only the gauge degrees of freedom, so we may 
set it to zero. 

We next calculate 

A3 = -E;^(0|«(^^)W|0) 

The connected contributions are given by 

A3 = -Ei(o|[k(^3n,«1^3lo)-E^(o|(^3)Mk^3],«1|o) 

-E^(0|k(^^)1 [Hha^m-^hQ\m)\a^] [a,//^]|0) (C.8) 

To calculate the terms on the second line, note that since depends on a and only 
in the combination (a + a^), we have 

[Hs,a'']^[a,H,] = J20E, 

where we have expanded the left hand expressions into a set of terms with distinct H2 
eigenvalues. Then 
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Using these expressions, we have 



^3 = -E;^(0|^k+i..^i^.+^J0)-E;^(0|(^k-i^.C?i^,-i^J0) 



where = | is the energy of the oscillator under consideration. 

To proceed, we note first that all disconnected terms in the expression above vanish 
since these will have = 0. To calculate the connected terms, we may use the explicit 
expressions for the operators O obtained from 



Note that the fermions terms contain only pieces with £'0 = and therefore will not 
contribute. The a, /3, and x contributions may be calculated separately using the 
propagators given above. For example, the (3 term will be 

2i 



-2 4f 



7V(7V2_i) (1)2 _ (1)2^6^ 

N{m - 1) J 

where in the second line the third factor is from the two ways of contracting the /?s, 
the fourth factor is from the propagators, and the final factor is from the sum over m. 
Similarly, the a contribution gives 

9 ^2(j + 2)2(2j + 3) 



N{N^ - 1) ij + 1 

Together, these sums may be evaluated to give 



N^{N + 1] 
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The contribution from the x terms is 



10^ 



N-l 



XXY 

3 7V(7V2 



A?^^ = E jU + 1) - -36 

^) j=0 



The remaining contribution to A3 comes from the terms on the first hne of equation 
(C.8), which we denote by A3. To compute these, we use 



[[a, Hs] , at] = 3J ^ ^^_ (a + at) (C.9) 



iV(iV2 _ 1) 

Then the first hne of equation (C.8) becomes 



= -3-3^ ^^J_^) ((0kg3|0) + (0|gV|0)) 

108 , 1 , 1 ^ 1 1 , 
-A - -B - -C + -D) 



iV(iV2-l)' 6 3 3 2 
27(2iV + l) 
m{N + l) 

Combining all the terms in A, we see that A^^^ cancels Af^-^^, while A|^^^^, A^^^ , 
and A3 sum to zero, so that 

A = 

Wc now proceed to calculate A{n = 2}. In this case, the H,^ contribution involves 
only the terms with each Y contracted with and a or at. Explicitly, we find 

A4 = ^(aaTV(-^[y^y^■][y^y^■])atat|o). 



1 1 

'24 
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N{N^ - 1) 



24Tr {[J\ J^Y) 



N{N^ - 1) 



where the ( • )* indicates that we are keeping only terms with no contractions between 
initial and final as. To compute the H^H^ terms, we note that 

2 o 
= «oo + ■ ■ ■ 



3^N{N^ - 1) 



where the additional terms have no more than one power of ctoo and therefore cannot 
contribute to A. Then 

A3 = n-^mHiym,/;)* 

= - E(qQ Q'oo Q'oo Q^Q^) _ 1) ^ 
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29N{m - 1) 
-270 



1 4 



(Jfj ^ > (3 • [36] - 1 • [36] - 3 • [144]) 



N{m - 1) 



where in the last hne, terms in {} brackets are from propagators, and the terms in [] 
brackets are combinatorical factors for the contributions with Ej, — —1/3, Ej, — 1/3 , 

and Ei = 1 respectively. 

Thus, for n = 1 the total energy shift is zero and for n = 2 we find 



We may now use (C.7) to write the result for the leading energy shift of the breathing 
mode with occupation number n, 



Thus, for the low-energy states, the leading perturbative correction is proportional to 
the bare coupling — {R/ (fiN)}^, which remains finite in the large N limit. 

D Height of the potential 

In this appendix, we will demonstrate explicitly the existence of a path between the 
irreducible vacuum (labeled by + 1) and the vacuum labeled by (A^, 1) along which 
the maximum energy is of order (/xiV/i?)^-^ = rather than 

Our choice of path is physically motivated by the following picture. The initial 
vacuum state corresponds to a single membrane sphere of radius while the final 
state corresponds to a slightly smaller sphere of radius ^^^^^^ plus a small sphere of 
membrane at the origin. If these were really continuum membranes, one way to move 
between these two points continuously at a seemingly small cost in energy would be to 
form a thin tube of membrane from the south pole of the original sphere to the origin 
where the small second sphere could grow from the tip of the tube. In this way, some 
of the membrane from the outer sphere could be transferred to the origin, yielding the 
desired final configuration, as shown in figure 2. 

We can translate this continuum picture to a path between the desired initial and 
final matrix configurations as follows. First, write down a map from the unit sphere 

^°In terms of the natural metric Tr {{X^)"^) on the space of matrices, this choice for the final state 
represents the closest vacuum. 



A = 



A3 + A4 
216 



N{m - 1) ■ 



A 
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(the membrane worldvolume) to the target space which is a function of some parameter 
e which will be related to the radius of the inner sphere. At a given value of e there will 
be some line of latitude near the south pole of the worldvolume sphere across which 
there is a singularity corresponding to the infinitesimal tube. The area above this line 
will map to the outer sphere, while the (small) area below this line will map to the 
small inner sphere. 

Once we have defined such a continuum mapping, we may obtain a corresponding 
path in the space of matrices by expanding the continuum map in terms of spherical 
harmonics, 

xHe,^/j,e)=J2xlJe)Yije,<l>) 
and then replacing the spherical harmonics with matrix spherical harmonics with I < 



^NxN 
■ Im 



KN 



Without explicitly specifying a map, it is clear that a natural choice is to preserve 
axial symmetry at all points along the path, so that we may choose to be a function 
of 9 only and + ix"^ = r(6')e*'^. In this case, the expansion of x^{9,(p) will involve 
only spherical harmonics with m — 0, while the expansion of x^ + ix"^ will involve 
only spherical harmonics with m — 1. As a result, at every point along the path, our 
matrices will take the form 



3R 



I ai 



V 



02 



I 



V 











since the matrices YJo are all diagonal, while the matrices Yw have nonzero entries only 
in the first row above the diagonal in each column. 

It is convenient to write the potential explicitly in terms of these variables. We find 

^ = (^)' (^'^"^ + - + E^/ - «^ + 1)') (D-i) 

Since this expression is quadratic in the variables 6^, it is not hard to determine the 
set of 6's for which the potential is minimum for a given set of a's. The extremum 
condition for h] is 



26/ 



-7 + ^?-i(^m 







Thus, at the minimum of the potential for a fixed set of a's, we will have for each I 
either 6;=0 or 

h]^^ - 2h] + hU = 2(a,+i -ai + lf-^ (D.2) 



^-'^Wc will not need to explicitly specify the mapping, but a natural choice would be to take the 
sum of the radii of the target space spheres to be a constant, and impose the condition that the map 
preserves ratios of areas. 
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Our problem is now reduced to specifying a/(e) such that the initial configuration 



IS 



N 



N 



--) 

2 ' 



and the final configuration is 



N -I 



,0) 



In principle, one could let the value of cat parameterize the path and explicitly extremize 
the potential in terms of the other a^'s to determine the lowest energy path satisfying 
our ansatz of axial symmetry. However, since the values of a/ for / < N change very 
little between their initial and final values it seems likely that this minimum path will 
have a maximum energy that does not differ significantly from that of the linear path 
between initial and final points. In other words, we expect that working to linear order 
in e will be sufficient to determine the order of magnitude of the barrier height. 
Thus, we take 



T + ^-^-2 



At e = 0, the values of h correspond to X+ — -^J'^ , 

k{N + l-k) 



N N 



so near e = wc will be on the branch with all &i's nonzero. Inserting the values of a 
into (D.2), we find 



/ -2 
1 

V 



1 



-2 / 



\ 



hi 



( 



I 



\ 



-2 



V -2 + |(A^ + 1) V + (A^ + l)e / 



The matrix M on the left may be inverted to give a symmetric matrix with 

iV + 1 -/ 



-k- 



N ■ 



k < I 



so we find that near e = 



bl 



k{N+l-k-e 



Note that 6jv in this formula is negative for e > {V2N + 3 — 1)/(A^ + 1), so for these 
values of e, we are no longer on the branch with all 6's non-zero. In fact, at e = 
{V2N + 3 - l)/(iV + 1) the values of 6, 



bl = k{N - k) 
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are already equal to their e = 1 values corresponding to J"*" in the (A^, 1) representation. 
It is not hard to check that these values, together with the values of a defined above, 
satisfy the equations (D.2) for all values of e. Thus, we may conclude that 



k{N + l-k-e-\{N+l)e^) e < 
' ^ k{N-k) e>VlMp 



Plugging the values for a and h into (D.l) we conclude that the energy along our path 
is given by 



V = 



,, 1 (Af+l) /i \2 , ^ V2N+3-1 

^ V -R / 648 Ar2 ^ ^ N+1 



Prom this expression, we find that the maximum of the potential along the path we 
have determined is 



R J [ 6487V V 

This represents an upper bound on the classical energy required to move between the 
vacuum labeled by {N+1} and the vacuum labeled by {N, 1}. It is clear that the same 
bound applies for a path between vacua {N + 1, A^i, . . . , A^„} and {1, N, Ni, . . . , Nn}. 
Thus, the minimum energy required to move from a given vacuum {A^i, A^2, ■ ■ ■ , N^} 
to another vacuum is not more than -rj^(u,N^in/R)^, where A^mm corresponds to the 
smallest sphere in the collection. 

Pinally, we note that by a series of these steps in which we move from a vacuum 
{iVi, • • • , Nn} to a vacuum {1, A^i — 1, N2, . . . , A^^}, and then by the reverse process 
to the vacuum {A^i, A^2 + 1, ■ ■ ■ , Nn}, we may transfer any number of momentum units 
between any two membrane spheres and in this way move between any two vacua of 
the theory. The maximum energy required for a given initial and final vacuum will be 
of order 

1 ffiN„ ^ 3 



E 



N \ B 

^ ' max ^ 



' max 



where Nmax corresponds to the largest sphere whose momentum is changed in the 
process. 
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